Elliptic flow in the Gaussian model of eccentricity fluctuations 
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We discuss a specific model of elliptic flow fluctuations due to Gaussian fluctuations in the ini- 
tial spatial X and y eccentricity components {((c"y — o"^)/(cr^ + o"y)) , {'2(Jxy / (o-'^ + ''"y))}- We flnd 
that in this model i'2{4}, elliptic flow determined from 4-particle cumulants, exactly equals the 
average flow value in the reaction plane coordinate system, (vrp), the relation which, in an approx- 
imate form, was found earlier by Bhalerao and Ollitrault in a more general analysis, but under the 
same assumption that V2 is proportional to the initial system eccentricity. We further show that 
in the Gaussian model all higher order cumulants are equal to W2{4}. Analysis of the distribution 
in the magnitude of the flow vector, the Q— distribution, reveals that it is totally defined by two 
parameters, «2{2}, the flow from 2-particle cumulants, and i'2{4}, thus providing equivalent infor- 
mation compared to the method of cumulants. The flow obtained from the distribution is again 

V2{4:} = (vrp). 

PACS numbers: 25.75.Ld, 25.75.-q 



I. INTRODUCTION 

Elliptic flow is an important observable in heavy ion collision experiments, which provides valuable information 
about the physics of the system evolution starting from very early times. Large elliptic flow values observed recently 



in experiments at RHIC 



;yst 
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are often used as an evidence for early system thermalization and as an argument for 
the creation of a new form of matter, sQGP, the strongly interacting quark-gluon plasma. With high statistics data 
obtained in the last few years at RHIC the analysis of elliptic flow becomes dominated by systematic uncertainties, 
mostly by inability to separate the so-called non-flow correlations (azimuthal correlations not related to the orientation 
of the reaction plane) and the effects of flow fluctuations Flow fluctuations can be due to different reasons: one 
that has attracted much attention recently is the fluctuations in initial eccentricity of the participant zone. Below we 



discuss only the flow fluctuations related to eccentricity fluctuations ^,3,151. In this paper we review the definitions 
of the different coordinate systems relevant to flow analysis. Then we discuss a particular model of eccentricity 
fluctuations. Within this model we show that by studying azimuthal correlations of produced particles at midrapidity 
it is in principle impossible to separate non-flow correlations from flow fluctuations effects as all observables contain 
the same combination of the two effects. 



* E-mail address: voloshin@wayne.edu (S.A. Voloshin) 
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FIG. 1; The definitions of the RP and PP coordinate systems. 



FIG. 2: The definition of the EP coordinate system. 



II. FLOW COORDINATE SYSTEMS 

We call the coordinate system defined by the impact parameter and the beam direction the reaction plane coordinate 
system, and use subscript RP to denote quantities in this system (see Fig. [T]). Then the orientation (azimuth) of the 
impact parameter vector in the laboratory frame is given by "^np. The principal axes of the participant zone will 
define the participant plane coordinate system with the corresponding angle ^'pp, and with the xpp axis pointing 
in the direction of the semi-minor axis of the participant zone. We use PP subscript for quantities defined in this 
system. 

The orientation of the flow vector Q = {Qx,Qij} = {X^i cos 20^, sin20i}, where the sum runs over all particles 
in some momentum window, defines the second harmonic event plane (see Fig. [2]) with corresponding azimuth 'I'bp, 
Qx = Q 0082^^ EP, Qy = Q sm2'^ EP- Although we use Q in this paper, in practice one would use q = Q/\/N in order 
to minimize the effect of the multiplicity spread within a centrality bin Q . For a given orientation of the participant 
plane. pp, anisotropic flow develops along this participant plane. 

The orientation of the participant plane can be also characterized by the eccentricity vector with coordinates 




part 



2(7x1, 



(1) 



y ' part 



where ct^ — {x'^) — {xY 
the participants in a given event 
magnitude, Spart = 




= {xy) — (y) (x) , and the average is taken over the coordinates of 
. The eccentricity vector direction is given by ^'pp = atan2(£j,, e^), and its 
epp, is called the participant eccentricity (see Figs.[3j|4]) in contrast with the reaction 



plane (or standard) eccentricity Ex = srp with its mean value defined to be 



{£rp) 



(2) 



This mean value is approximately Sopt, the optical eccentricity determined by the optical Glauber model 
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III. GAUSSIAN MODEL FOR ECCENTRICITY FLUCTUATIONS 



In events with fixed e, botli in magnitude and orientation, the flow vector on average points along e, but with the 
magnitude and orientation of the flow vector fluctuating due to finite multiplicity of particles used in its definition. 
As can be seen from simulations using the MC Glauber model ^, Q| in Fig. [5l the distributions in Sx and Sy are 
well approximated by a Gaussian form with widths approximately equal in the two directions. There exists some 
deviation from a Gaussian form in peripheral collisions, but even there the deviations are small, so we proceed with 
the Gaussian ansatz. We denote the equal widths in Sx and Sy by cr^. The distribution in the magnitude of the 
eccentricity, Spart, can be obtained by integration over angle of the vector e as a two-dimensional Gaussian (see, for 
example, the derivation in Q), and is given by 

dn e part J f £part{£B^p)\ ( 4<"t + (^flp) ^ \ I \ \ I ■3\ 

^ = [ ,2 } y J - BG(s,.,.; {snp) , ex.), (3) 

where we have introduced a short hand notation BG(a;; x, a) for the "Bessel-Gaussian" distribution with one variable 
argument and two constant parameters (see Fig. [6]). Note that in BG(£part; {^rp) , <^e): £part is an eccentricity as given 
in PP but (epp) and describe the 2-D Gaussian distribution in the RP system. The distribution is normalized to 
unity. For later use we provide a few moments of the distribution BG(a;; x, a), where a; is a generic variable (not the 
X-axis) : 



2a V 4cr2 / V 2 



(4) 



(x^) = x^ + 2a', (5) 

(^4) = x* + 8x^a^ + 8a\ (6) 

(x^) = x^-f 18x^2 72x^(74 + 48cr6. (7) 

Note that the parameter a is not the variance of this distribution; the latter would be given by 

al = {x')~{x)\ (8) 

with (x^) and (x) given above. Also, from Eqs. ^ and ^ it can be shown that 

2(x2)'-(x4) =x4 (9) 
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FIG. 5: (top) Distribution in together witli Gaussian fits for (left to right) central, mid-central, and peripheral collisions, 
(bottom) The y and x distributions for the mid-central case. All curves have been normalized to the same area. 

and 

{x'')-<i{x^){f) + \2{ff ^\x\ (10) 

In very central collision, the non-zero eccentricity of the overlap region is defined mostly by fluctuations, {epart) 3> 
(erp). This limit corresponds to a; <C cr in Eqs. ([3]-[5]). One finds in this limit (x) — a^fiTpi and Oxj (x) — ^/A/n ~ 1, 
the relation first derived in [sj. 

Figure [6] shows the distribution in Epart from the MC Glauber calculation, together with the fit to the BG form. 
The quality of the fit is good, and the extracted fit parameters shown in Table U agree well with those extracted 
directly from the distributions of Fig. [5] bottom for Ex and Ey. 
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FIG. 6: Distribution in Spart for mid-central collisions (4 < 6 < 6 fm) and fit to the BG shape. 

TABLE I: Comparison of a Gaussian distribution of e in the RP system with the Bessel-Gaussian fit in the PP system for 
mid-central collision. 

G, e^, (Fig. El) 0.1384 ± 0.0001 0.0935 ± 0.0001 
G, Ey, (Fig. [5} 0.0000 ± 0.0001 0.0923 ± 0.0001 
BG (Fig. El) 0.1344 ± 0.0002 0.0957 ± 0.0001 



IV. FLOW FLUCTUATIONS IN A GAUSSIAN MODEL OF ECCENTRICITY FLUCTUATIONS 

We start our consideration by deriving the flow vector distribution. One can approach this problem starting from 
two different coordinate systems: the participant coordinate system or the reaction plane one (see Fig. mi. In the PP- 
system the y coordinate of the flow vector is not affected by flow (and/or flow fluctuations), only the x component 
is, which might be taken as a simplification. On the other hand the fluctuations in participant eccentricity (and 
correspondingly, in flow) have the BG form, which is more difficult to take into account analytically. Somewhat easier 
(though, obviously, equivalent) is to perform the analysis of the Q-distribution in the reaction plane system. In the 
i?P-system both components of the flow vector are affected by eccentricity fluctuations, but the fluctuations arc of 
Gaussian form, with the same widths in the x and y directions. Assume that on average, flow is proportional to 
eccentricity with proportionality coefficient k: 

V2 = KEparf (H) 

For events with fixed e = {ex,ey} this leads to {Qx)^ ^ NkEx, {Qy)^ = NnSy. For the overall distribution one finds 
that the fiow vector is a two-dimensional Gaussian distribution with {Qx) = Nk {srp) , (Qy) = 0, and widths in the 



two directions given (see Q], [l^) by 

'^Qy^^^Es"^^^,^ ^ = \m^-{cosm,)) + {N~l){2K'a',+6)], (12) 

where N is the number of particles, and S is the non-flow contribution defined by (mm*) = {cos{2(pi — 2(j)j)) = u| + (5, 
with M being the single-particle unit (second harmonic) flow vector. Neglecting the contributions of the fourth 
harmonic flow and the {k {£x)Y term, both less than or of the order of 10^^ - 10~* compared to unity, one finds that 
the widths in both directions are the same: 

4. = ^Qy = ^^[1 + - 1)(2«^''^? + 5)], (14) 

Note that k {£b,p) ~ {urp) = v gives the real flow as calculated with respect to the reaction plane and the standard 
deviation of v along the reaction plane axis is na^ = a^x- The distribution in flow vector magnitude would be given 
then by 

dn/dQ^BG{Q;NK{£Rp) .UQx). (15) 

Let us now calculate from 2-particle and four-particle cumulants , M2{2} and M2{4}, using the Gaussian 

ansatz for flow fluctuations. 

M2{2}' = (cos(2(/., - 20,)) = {vf) + <5 = Ac^ (4^,,) + 6. (16) 

Using Eq. ([5]) this becomes 

Similarly, for the fourth order cumulant result, using Eq. ([9]), 

M2{4}^ = 2 (cos(2(/), - 24>j)f - (cos(20, + 20, - 2<t>k - 20,„)) = 2 («|)' - (m|) = = (mj^p)^ . (18) 

Note that in this approach (Gaussian ansatz) W2{4}^ is always well defined as the cumulant does not change sign. In 
our model the relation (|18p is exact, but in an approximate for m ( and using a different treatment of the eccentricity 
fluctuations) it was derived earlier by Bhalerao and OUitrault [12|, who were the first to note that the fourth order 
cumulant flow measurements are mostly unaffected not only by non-flow effects but also by flow fluctuations. 
Proceeding further, for the difference of the two cumulant results one obtains from Eqs. (|17|1 and ([T8|) 



V2{2Y - V2{^Y - 2^2 (jf + 5 = 2al^ + S, (19) 



unfortunately the same parameter that defines the Q distribution width in Eq. ()14p . The last observation rules out 
(in the Gaussian ansatz) the possibility to measure both fluctuations and non-flow by combining information from 
Q-distributions and cumulants. Neither do higher order cumulants provide new information. Using Eq. (jlOp one flnds 
out that 

M2{6}6 = ((^2') - 9 (vt) {vl) + 12 {vlf) /4 = {vnpf . (20) 
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One can show that in this model all higher order cumulants are given by the corresponding power of (vrp). Another 
way to look at this is to apply Eqs. ([9]) and (fTO|) directly to the Q distribution Eq. (fTSjl . One finds that the 
combinations usually associated with flow cumulants are given by corresponding powers of Nvrp, for example 
2{Q^)-{Q^) = {NvnPr. 



V. FITTING Q-DISTRIBUTIONS 



As can be seen by comparing Eqs. ([T4|) and (fT9|l . V2{2} and U2{4} completely define the form of the Q-distribution, 
and can be used as an alternative set of parameters compared to that in Eq. (|15[) . If one tries to fit the Q-distribution 
with a functional form determined by three parameters, e.g. (v), ay, and S, these parameters should satisfy the values 
of W2{2} and W2{4} (which provides only two equations), and all three can not be determined. 

There can be different functional forms used to describe flow fluctuations along the PP axis. Most often used 
are the Gaussian form G(u; (v) ,(7y) and the Bessel-Gaussian BG(w;wo, cr) discussed above. Both of them have two 
parameters, which, as we know can not be determined separately, so they must be correlated. 

Assuming the JiG{v;vo,a) form for flow fluctuations to fit the Q-distribution, which would correspond to a two- 
dimensional Gaussian distribution in the reaction plane coordinate system, one would find from Eqs. (|17p and (|18p 
that the parameters are correlated according to 

V2{2}'^ = const = + 2a^ + S (21) 
U2{4} = const = vq. (22) 

The mean and the variance of the v distribution would be given by Eqs. ^ and ([S]), but since a can not be determined 
independent of 5, (v) is also undetermined. 

If one uses the Gaussian form for flow fluctuations in the PP-system, one would find that the parameters are 
correlated according to 

V2{2}^ = const =^ {vf + al + S (23) 
«2{4}^ = const = ^{vf-2{vfa^,-at « {vf ~ (24) 

or equivalently 

V2{2f - V2{Af = const K 2al + 5. (25) 

The above two equations are derived in the approximation of CTu ^ (u) but for Gaussian fluctuations in v the exact 
formula can be used. Again, as a can not be determined independently, (v) is also undetermined. 



VI. SUMMARY 



We find that in the Gaussian ansatz, fitting Q-distributions does not bring any more information than that provided 
by cumulants. It is not surprising - if the distribution is defined just by two parameters one can not get more than 
?;2{2} and W2{4} already provided. Note that under this ansatz all the higher order cumulant v values are the same. 
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The origin of the "problem" can be traced to the Gaussian ansatz. It is known that for a Gaussian distribution all 
the cumulants higher than rank two arc zero. The latter means that if the collective fluctuations are of the Gaussian 
type one can never prove that the fluctuations exist by any type of correlation analysis using only particles under 
consideration (no external information) . A similar problem was observed earlier in a temperature fluctuation study of 
many-particle transverse momentum correlations IJ]. Unfortunately, deviations from a Gaussian distribution might 



be too small to observe. Such deviations would show up in the bad quality of the Q-distribution fits based on the 
Gaussian ansatz, or in a small differences between higher order cumulant v values. 

The fact that all higher order cumulants are the same and determined by the value of flow in the reaction plane (not 
the participant plane), and that that fitting of Q-distribution yields the same value, explains the consistency between 
t;2{4} and t-2{ZDCSMD} Q, 

which is calculated with ZDC-SMD as event plane and is supposed to be sensitive to 
V2 in the reaction plane, as well as the consistency between f2{4} and V2{Q — dist} 

Ref. H used a model of fiow fiuctuations in which fiow fiuctuates only in the impact parameter direction. The use 
of the detectors which measure spectator neutrons, advocated in 1^, is justified for that model, but would yield zero 
results for the case of fluctuations discussed in this paper. 

The authors thank P. Sorensen for fruitful discussion. This work was supported in part by the HENP Divisions of 
the Office of Science of the U.S. Department of Energy. 
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